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ABSTRACT 

We detect a Neptune mass-ratio (5 ~ 8 x 10~^) planetary companion to 
the lens star in the extremely high-magnification (A ~ 800) microlensing event 
OGLE-2005-BLG-169. If the parent is a main-sequence star, it has mass 
M ~ 0.5 implying a planet mass of ~ 13 and projected separation 
of ~ 2.7 AU. When intensely monitored over their peak, high-magnification 
events similar to OGLE-2005-BLG-169 have nearly complete sensitivity to Nep- 
tune mass-ratio planets with projected separations of 0.6 to 1.6 Einstein radii, 
corresponding to 1.6-4.3 AU in the present case. Only two other such events 
were monitored well enough to detect Neptunes, and so this detection by it- 
self suggests that Neptune mass-ratio planets are common. Moreover, another 
Neptune was recently discovered at a similar distance from its parent star in a 
low-magnification event, which are more common but are individually much less 
sensitive to planets. Combining the two detections yields 90% upper and lower 
frequency limits / = O.ST^q'I^ over just 0.4 decades of planet-star separation. In 
particular, / > 16% at 90% confidence. The parent star hosts no Jupiter-mass 
companions with projected separations within a factor 5 of that of the detected 
planet. The lens-source relative proper motion is ~ 7-10 mas yr~^, implying 
that if the lens is sufficiently bright, / < 23.8, it will be detectable by the Hubble 
Space Telescope by 3 years after peak. This would permit a more precise estimate 
of the lens mass and distance, and so the mass and projected separation of the 
planet. Analogs of OGLE-2005-BLG-169Lb orbiting nearby stars would be dif- 
ficult to detect by other methods of planet detection, including radial velocities, 
transits, or astrometry. 

Subject headings: gravitational lensing - planetary systems - Galaxy: bulge 

1. Introduction 

The regions of our Solar System beyond Mars contain giant planets, as well as asteroids 
and comets, which are believed to be remnants of the process that formed these bodies. 
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and also moons, which may serve as analogs to bodies involved in the late stages of planet 
formation. However, despite some 170 planet discoveries over the past decade, the analogous 
regions around other mature stars remain relatively inaccessible to us. As radial velocity 
(RV) survey time baselines have grown, they have begun to detect gas giants in these regions, 
but RV is sensitive to Neptune-mass planets only when they are much closer to their parent 
stars. Transit surveys are even more heavily biased toward close-in planets. Astrometric 
sensitivity does peak at large orbits but is fundamentally restricted to orbital periods that 
are shorter than the survey. 

By contrast, microlcnsing sensitivity peaks at the Einstein ring, which is typically at 
2-4 AU, depending on the mass and distance of the host star, and it extends several times 
farther out. Moreover microlensing detections are "instantaneous snapshots" of the system, 
not requiring an orbital period to elapse. Hence, microlensing can potentially yield important 
information about the intermediate-to-outer regions of extrasolar planetary systems that are 
difficult to probe by other techniques. 

Gravitational microlensing occurs when a "lens" star becomes closely aligned, within 
an angular Einstein radius ^e, with a more distant source star. The source is magnified by 
an amount that grows monotonically as it approaches the lens, and diverges inversely with 
separation for extremely close encounters (Einstein 1936; Paczynski 1986). 

Planets hosted by the lens star induce two, generally distinct, perturbations on the 
single-lens magnification pattern: a small "planetary caustic" (closed contour of formally 
infinite magnification) directly associated with the planet and an even smaller "central caus- 
tic" closely aligned with the host. All microlcnsing events have the potential to probe the 
planetary caustic, but owing to its small extent and the random trajectory of the source, the 
probability for the source to encounter the perturbed region is low. 

By contrast, the small fraction of events that reach very high magnification (and so 
small lens-source separation) automatically probe the region of the central caustic (Griest 
& Safizadeh 1998; Rhie et al. 2000; Bond et al. 2002; Rattenbury et al. 2002). Hence, even 
low-mass planets lying anywhere sufficiently near the Einstein ring are virtually guaranteed 
to perturb the light curve (Abe et al. 2004). Since the peak of the event can, at least in 
principle, be predicted in advance, it is possible to focus limited observing resources to make 
the intensive observations over the peak that are required to detect and characterize the 
relatively subtle signal. 

The Microlensing Follow Up Network^ (/^FUN) has adopted a strategy of trying to rec- 
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ognize the few prospective high-magnification events among the roughly 500 microlensing 
events annually alerted by the Optical Gravitational Lens Experiment^ (OGLE-III) Early 
Warning System (Udalski 2003) and the roughly 50 events annually alerted by the Microlens- 
ing Observations for Astrophysics^ collaboration (MOA). During the 2005 season, this strat- 
egy led to the discovery of 2 planets. The first was the Jovian mass-ratio OGLE-2005-BLG- 
071Lb (Udalski et al. 2005). Here we report the second such planet, a Neptune-mass-ratio 
companion to the lens in OGLE-2005-BLG-169. 

2. Event Discovery and Data 

On 2005 Apr 21, the OGLE collaboration alerted OGLE-2005-BLG-169 as probable 
microlensing of a faint (/ = 19.4) source toward the Galactic bulge, using the 1.3m Warsaw 
telescope in Chile (operated by the Carnegie Institution of Washington). After observa- 
tions by OGLE and the 1.3m //FUN SMARTS telescope in Chile showed the event to be 
extremely high magnification, the observers at the 2.4m //FUN MDM telescope in Arizona 
interrupted their regular program to obtain more than 1000 exposures over the peak. Addi- 
tional data come from the 0.35m /xFUN Nustrini telescope in Auckland, New Zealand and 
the 2.0m PLANET/RoboNet^ Faulkes Telescope North in Hawaii. We analyze a total of 
(340,22,1025,74,31) images with typical exposure times of (120,300,10,120,100) seconds in 
the (/, /, /, clear, R) passbands, respectively, from these 5 telescopes. In addition, 6 V band 
images from //FUN SMARTS permit determination of the source color. 

All data were reduced using the OGLE data pipeline based on difference imaging analysis 
(DIA) (Wozniak 2000). The /iFUN MDM data were also reduced using the ISIS pipeline 
(Alard & Lupton 1998; Alard 2000; Hartman et al. 2004). To test for any systematics in 
this crucial data set, we report results below derived from these two completely independent 
pipehnes. See also Figure 1. 

3. Light Curve Model 

Planets are discovered in microlensing events from the brief perturbation they induce on 
a single-lens hght curve (Mao & Paczyhski 1991; Gould & Loeb 1992): F{t) = FsA[u{t)]+Fy,, 



^http://www.astrouw.edu.pl/~ ogle/ogleS/ews/ews.html 
•^http: / / www.massey.ac.nz/~iabond/alert / alert.html 
*http://planet. iap.fr/, http://www.astro.livjm.ac.uk/RoboNet 
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where A{u) = {u'^ + 2)/u{u^ + 4)^2^ u{t) = (r^ + r = - io)AE, Here, F(i) is the 

observed flux, is the source flux, Fb is flux due to any unlensed background hght, A is 
the magniflcation, and u = (r, Mq) is the vector position of the source normahzed to (6'e), 
expressed in terms of the three geometrical parameters of the event: the time of closest 
approach to, the normalized impact parameter mq, and the Einstein timescale ^e- 

To describe the planetary perturbation, 3 "binary-lens" parameters are required in ad- 
dition to the 3 single-lens geometric parameters (toj Mq, ^e)- These are the binary mass ratio, 
the separation of the components (in units of ^e), and angle of the source trajectory 
relative to the binary axis, a. Finally, a seventh parameter, p = 9^/6e, is required whenever 
the angular radius of the source 6** plays a significant role, in particular, whenever the source 
crosses a caustic. 

For planetary lenses with caustic crossings, there are 7 directly observable and pro- 
nounced light-curve features that directly constrain the 7 model parameters, up to a well- 
understood two-fold degeneracy (Dominik 1999) that takes b <-> b~^. Three of these, the 
epoch, duration, and height of the primary lensing event, strongly constrain {to, tE, Uq). The 
remaining four, the two caustic-crossing times (entrance and exit) and the height and dura- 
tion of one of the caustic crossings, then constrain (6, q, a, p). 

As we now show, OGLE-2005-BLG-169 is indeed a caustic-crossing event. However, 
only the caustic exit (but not entrance) was well-resolved. This leads to a 1-dimensional 
degeneracy among the model parameters, which is then partially broken by secondary, less 
pronounced, features of the light curve. 

The residuals to the single-lens models shown in Figure 2 exhibit a kink in slope at 
At = 0.092 day, where At is the time elapsed since HJD 2453491.875 (2005 May 1 09:00). 
The kink is equally apparent in both sets of MDM reductions. While the incident slope 
of this kink in the light curve depends on the particular single-lens model used, the slope 
discontinuity is model-independent. Such a change in slope is induced by a caustic exit, 
when the trailing limb of the source crosses the caustic, causing two of its images to merge 
and finally vanish. 

To extract model parameters, we undertake a brute-force search of parameter space. 
We hold the three parameters (h, a, q) fixed at a grid of values while minimizing over the 
remaining four parameters, using the values of (to,'Wo,tE) derived from the overall shape of 
the light curve as seeds. We find two distinct local minima that obey the b ^ b^^ degeneracy. 
These minima are embedded in elongated Ax^ valleys (Fig. 3), with [a (in radians): b] axis 
ratio ~ 100, which occur because there are only 6 pronounced features in the light curve to 
constrain 7 model parameters. However, while the caustic entrance is not well resolved, the 
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point at At — —0.1427 day does lie on this entrance and thereby singles out the solutions 
shown in Table 1. Nevertheless, to be conservative, we quote 3 a errors for light-curve 
parameters, which take account of both the elongated valleys and the two reductions. Most 
importantly, 

q = 8tlx 10-^ b = 1.00 ± 0.02 (3 a). (1) 



4. Mass of the Host Star and Planet 

In microlensing events, the physical parameters (lens mass M, lens-source relative par- 
allax TTrei, and relative proper motion /x) are related to the "observable" event parameters 
(^E, ^E, and the microlens parallax (Gould 2004), tte), by 



9E^^/KM^T^l, tte = \ , = — , (2) 

V TTrel 

where k. = 4G/(c^AU). Of these, only is routinely measurable. However, in caustic- 
crossing events, one can usually also measure 9e, which both directly constrains the mass- 
distance relation and yields a measurement of /i. We first determine the angular source radius 
9^ using the standard approach (Yoo et al. 2004a), finding 9^ — 0.44 ± 0.04 /xas. Together 
with the parameter measurements p — 4.4lQ;g x 10~^ and ^e = 43 ± 4 days, this yields 3 a 
ranges, 

^E = — = 1.00 ±0.22 mas, // = — = 8.4 ± 1.7masyr-\ (3) 

P tE 

Another constraint comes from the upper limit on the lens flux, which cannot exceed 
the background flux measurement, Fb (corresponding to Jb = 19.8). Wc also derive from the 
light curve a weak constraint on the microlens parallax, tteji = —0.086 ± 0.261, which we 
include for completeness. Here, 7rE,i| = ttecos^/' and ip is the angle between the direction of 
lens-source relative motion and the position of the Sun at to projected on the plane of the 
sky (Gould 2004). 

These measurements, together with a Bayesian prior for the lens distances and masses 
from a Han & Gould (2003) Galaxy model and Gould (2000) mass function (but with a 
Salpeter [—2.35] slope at the high end), together with the assumption that all stellar objects 
along the line of sight are equally likely to harbor planets, yield a probability distribution 
for the lens mass. We find that the probabihties that the lens is a main-sequence star (MS), 
white dwarf (WD), neutron star (NS), and black hole (BH) are respectively 55%, 32%, 11%, 
and 2% . However, Neptune mass-ratio planets around NSs must be quite rare because the 
~ 0.2 second oscillations they induce would easily show up in pulsar timing residuals. Since 
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BHs and NSs form by a similar process, this may also argue against BHs as the host. MSs 
and WDs represent different life stages of the same class of stars. If the host is a MS, then 
the median and 90% confidence interval for the mass and distance are given by 

M = O.49l°;|M0, = 2.7^^1 kpc, (90% confidence) (4) 

implying that the best estimates for the planet mass and separation are rUp — qM ~ 13 
and rx = bO^Di^ ~ 2.7 AU, where is the distance to the lens. Since WD masses are sharply 
peaked at M ~ 0.6 Mq, the corresponding planet characteristics lie near the center of the 
same range. Future observations by the Hubble Space Telescope could distinguish between 
MS and WD hosts, as well as measuring the mass and distance of the former provided that 
the lens is brighter than I ~ 23.8 (Bennett & Anderson 2006). 

5. Limits on other Companions 

Does lens OGLE-2005-BLG-169L have planets other than OGLE-2005-BLG-169Lb? 
Central-caustic events allow one to address this question because all the companions to 
the lens star perturb the central caustic (Gaudi et al. 1998). Moreover, the combined per- 
turbation from two planets is very nearly the sum of the separate perturbations, unless the 
two planets are closely ahgned (Rattenbury et al. 2002; Han 2005). This allows us to subtract 
the perturbation from the one detected planet and apply the same search technique to the 
resulting (nearly single-lens) light curve to look for others. We find none, and so place upper 
limits on the presence of other planets. In particular, we exclude Jupiter-mass (g = 2 x 10~^) 
planets at projected separations within a factor 5.5 of the Einstein radius (0.18 < b < 5.5) 
and Saturn- mass planets within a factor 3.5. 

6. Discussion 

6.1. Low-mass Planets Are Common at Several AU 

OGLE-2005-BLG-169 is one of only three high-magnification events with sensitivity 
to the central caustics induced by cold Neptunes. It and one other event, MOA-2003- 
BLG-32/OGLE-2003-BLG-219 (Abe et al. 2004), were sensitive to g = 8 x 10"^ planets 
throughout the "lensing zone" (0.6 < 6 < 1.6), while OGLE-2004-BLG-343 (Dong et al. 
2006) was sensitive over 40% of the lensing zone. This implies (/) = 1/2.4 = 42% for the 
expected fraction of stars hosting cold Neptunes within the lensing zone (0.4 decades of 
projected separation), but with very large uncertainty due to small number statistics. 



-9- 



However, we can improve our estimate by incorporating the detection of OGLE-2005- 
BLG-390Lb, which is a very similar {b = 1.6, q = 7.6 x 10~^) planet that was detected through 
the other (planetary-caustic) microlensing channel (Beaulieu et al. 2006). The expected 
number of detections of g = 8 x 10^^ lensing-zone planets (through this planetary-caustic 
channel) was 1.75/ for 1995-1999 (Gaudi et al. 2002). We estimate that the expected number 
for 2000-2005 is a factor 1.5 higher based on a comparison of the survey characteristics during 
these two periods. Hence, the single detection through this channel yields (/) = 23%, also 
with large errors. If we combine the two channels, and impose a uniform prior, we find a 
median and 90% upper and lower limits oi f — O.ST^oii) particular, a 90% confidence 
lower limit of / > 16%. 



6.2. Unique capabilities of microlensing 

Can analogs to OGLE-2005-BLG-169Lb (with the same mass ratio and semi-major axis) 
be detected by other techniques, such as RV, transits, and astrometry? To make a strict 
comparison, it would be necessary to know three additional parameters of the event, namely 
the lens and source distances, Di^ and Ds, and the projection angle ip of its semi- major axis 
a to the line of sight, i.e., sin'^ = r±/a where r± is the projected star-planet separation. The 
stellar mass and semi-major axis are then related to the event parameters by M = 6^/ ktTi-ci 
and a = Di^O^b esc ip. Hence, assuming circular orbits, the velocity amplitude vsini and 
orbital period P would be 



. . qcsini /^EsinV' „ _i / sin^' . . 

.sm. = ^y^^^ 0.85ms y^^^^^-n., (5) 

^ V63esc3^(l-x)^, ^ 9.4 yrs-^ f ^) ■ (6) 
c ^ v^v y y 2.7kpcV 0.66 V 1.3 / ^ ^ 

Here x = Di^/ Do, and the evaluations are for = Imas, q — %x 10~^ and 6 = 1. These 
equations imply that RV detection of OGLE-05-BLG-169Lb analogs would be extremely 
difficult both because of the low velocity amphtude and the long period. Similarly, the long 
period renders transit detection impossible unless a transit experiment of considerably longer 
duration than Kepler were organized. Finally, the amplitude of astrometric motion would 
be 

« = ^E&cscV^,^ 28/xas — — -, (7) 

-'-^cmalog KpC i/analog ^-^ 

where Danaiog is the distance to the local analog system. This amplitude would be easily 
detectable by the Space Interferometry Mission (SIM), although at the typical periods, this 
would require an extended SIM mission. See Figure 4. 
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6.3. Neptunes or Failed-Jupiters? 

Both OGLE-2005-BLG-169Lb and OGLE-2005-BLG-390Lb arc in the cold outer regions 
of their planetary systems, with expected surface temperatures of ~ 70 K and ~ 50 K, 
corresponding roughly to the environments of Saturn and Neptune, respectively. Both have 
Neptune/Sun mass ratios, with absolute mass estimates of 13 and 6M®, respectively. 
They must have a large fraction of rock and ice, but whether these are covered with a thick 
coat of gas, like Uranus and Neptune, or whether they are "naked" super-Earths such as are 
theorized to have formed the cores of Jupiter and Saturn, is unclear. If such cores formed 
routinely but usually failed to accrete the ambient gas before it dispersed, this would account 
for the high frequency of these objects (Ida & Lin 2005). Moreover, the absence of gas giants 
within a factor 5.5 of the OGLE-2005-BLG-169's Einstein ring is consistent with the idea 
that the detected planet is such a "failed Jupiter" (Laughlin et al. 2004). One could gain 
further clues by mapping out the mass and separation distributions of a larger sample. 
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(DD,AG,RP); Polish MEiN 2P03D02124, NSF AST-0204908, NASA grant NAG5-12212 
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Table 1. Light Curve Models 
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.05 .1 
At [HJD -2453491.875] 



Fig. 1. — Light curves of 6 reference stars (black) compared to microlensed source (red) as 
derived by the DIA reductions of the MDM data. The deviation of this curve from a single- 
lens microlensing fit (derived from the non-MDM points) is shown in green. Horizontal blue 
lines are shown to aid in judging the constancy of the various stars. 
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At [HJD - 2453491.875] 



Fig. 2. — Data and best-fit model of OGLE-2005-BLG-169. Lower panel shows difference 
between this model and a single-lens model with the same (tQ,uo,tE, p)- It displays the 
classical form of a caustic entrance/exit that is often seen in binary microlensing events, 
where the amplitudes and timescales are several orders of magnitude larger than seen here. 
MDM data trace the characteristic slope change at the caustic exit {At = 0.092) extremely 
well, while the entrance is tracked by a single point (At = —0.1427). The dashed line 
indicates the time to. Inset shows source path through the caustic geometry and indicates 
the source size, p. 
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Fig. 3. — Ax^ contours relative to the respective minima for light-curve fits using the DIA 
and ISIS reductions of the MDM data, a = a/A^^ < 1,2,3,4,5,6,7 are shown by black, 
red, yellow, green, cyan, blue, and magenta. Note that abscissae are compressed, so the 
6a contours have an axis ratio of approximately 100 (with a expressed in radians), which 
refiects the 1-D degeneracy discussed in the text. Both reductions have their minimum at 
(a, g) ~ (120°, 8 x 10~^). Similar contours for 6 < 1 yield additional solutions that are 
included in Table 1. The full 3 a mass-ratio range is confined to 5-10x10"^ even allowing 
for both reductions. 
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Fig. 4. — Exoplanet discovery potential and detections are shown as functions of planet 
mass and semimajor axis. Potential is shown for current ground-based RV (yellow) and, 
very approximately, microlensing (red) experiments as well as future space-based transit 
(cyan), astrometric (green), and microlensing (peach) missions. Planets discovered using the 
transit (blue), RV (black), and microlensing (magenta) techniques are shown as individual 
points, with OGLE-2005-BLG-169Lb displayed as an open symbol. Solar system planets are 
indicated by their initials for comparison. 



